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'■^ ' Abstract 

We analyze here different forms of fractional relaxation equations of order 
rf^ ' V £ {0, 1) and we derive their solutions both in analytical and in probabilistic 

forms. In particular we show that these solutions can be expressed as cross- 
ing probabilities of random boundaries by various types of stochastic processes, 
, ^ which are all related to the Brownian motion B. In the special case ly = 1/2, the 

Mm . fractional relaxation is proved to coincide with Pr {suPq<^<j B{s) < U}, for an 

Mh ' exponential boundary U. When we generalize the distributions of the random 

(— I ^ boundary, passing from the exponential to the Gamma density, we obtain more 

and more complicated fractional equations. 
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1 Introduction 
in 

CN ' The following differential equation 

!>■ d 

o : -p{t) = -xpit), t > (1.1) 

^^ at 

is known in the physics literature as the relaxation equation. The solution to (jl.l[) . 

with initial condition p{0) = 1, is clearly equal to p{t) ~ e~'^*. Since the end of the 
k>( \ Nineties an intensive research activity has been developed, aimed at the application of 

; I ■ fractional calculus to mathematical physics: many classical equations have been niod- 

C^ I ified by substituting the integer-order derivatives with the fractional ones. Equation 

p.ip has been extended in the following fractional sense: 

|J;^(i) - -AV'(t), t>0 (1.2) 

where i^ e (0, 1) and -^ represents the fractional derivative according to the Caputo 
definition, i.e. 

^* I l^^W' for 7^ = 771, 
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with 771 = [aj + 1. Obviously, for 1^ = 1 the fractional relaxation equation (|1.2p 
coincides with the standard equation (|l.ip . 

Equation (|1.2p has been studied in some papers, such as [H], [H] and its solution 
was given analytically in terms of the Mittag-Leffler function as: 

Mt) ^ E.^ii-Xt''), (1.4) 

where 

£;a,M^)-ET^r^Tl^' «,/3eC, i?e(a),i?e(/3)>0. (1.5) 

f-'o r(ar + 13) 

The analysis of the fractional relaxation equation has mainly physical motivations, for 
instance to study the electromagnetic properties of a wide range of materials (which 
display a long memory, instead of exponential, decay, see [5S] and [55]) as well as the 
rhcological models for the description of some viscoclastic materials (see [TH], |H], [5D] 
and i27jV 

Moreover, the so-called Mittag-Leffler distribution has been often applied to statis- 
tics (for example in [T3] and [23]) or to queuing theory in [55]. 

Actually the solution ip^(t),t > can be expressed in probabilistic terms in two 
interesting forms, that we will present and explore here. The first form represents 
the probability of no events up to time t (or survival probability), for the so-called 
fractional Poisson process Af^{t),t > (see, among the others, [12], [30], [15], [1], and 
[3]). Indeed the following equality holds 

^At)^Poit)^P^WAt)^0} (1.6) 

and thus we can apply to '>p,y{t) the results obtained in the above cited articles. For 
example we will resort to the equality of the one-dimensional distribution between AC 
and a composition of the standard Poisson process N{t) with a random time-process 
7I/(i), i.e. N{Tu{t)),t > 0. Thus, thanks to ()1.6p . we can write 



OO 



Mt)^ e-^yq,{y,t)dy = FY{%{t)<U}, (1.7) 

Jo 

where q^{y,t) is the density of 71/ (which is itself solution to a fractional diffusion 
equation) and U is an exponential random variable with parameter A > 0. Formula 
()1.7|) is particularly interesting in the special case where v = 1/2, since it becomes 

A/2{t)= e'^y^^^dy = PT{\B{t)\<U}, (1.8) 

Jo VTTi 

where iJ is a Brownian motion starting from zero and with variance 2t. 

As a consequence, a second probabilistic interpretation of the solution to the frac- 
tional relaxation equation (|1.2[) can be given in terms of crossing probability of a 
random boundary by a standard Brownian motion, for i/ ~ 1/2. Indeed it is well 
known that the following relationship holds; 



Pr{|S(i)| <z} = Pr<^ sup B{s) < z )■ =Pr{B(.s) < z, V s G (0,i)}, 

Lo<,s<t J 

where the last expression is commonly referred to as crossing probability. 

For other values of v, an analogue result holds true, but for less known processes, 
such as the iterated Brownian motion (for v = 1/2") or the Aii^y process (for v = 1/3). 

Moreover, the expression (|1.7p shows that the solution to (|1.2p can be expressed 
as a standard relaxation with random time represented by Tu, i.e. as 'ip{T^{t)). The 



results given in |19| permit also to express the solution as a time-changed relaxation 
via an inverse stable subordinator E{t), i.e. as ipu{t) = ^{E{t)). In fact tlj{Tv{t)) and 
tj}{E{t)) share the one-dimensional distributions and therefore the two approaches can 
be considered equivalent. 

In the successive sections we analyze some extensions of the result (|1.7p in the 
following directions: 

• Wc consider other random time-processes in place of 71/ and therefore in p. 81) 
instead of Brownian motion: for example, the sojourn time of a Brownian motion 
on the positive half-line, the first-passage time of a Brownian motion through a 
certain level, the elastic Brownian motion (by analogy with the analysis carried 
out for the fractional Poisson process in [S] ) . 

• We consider a different random variable (i.e. the Gamma) instead of U in (|1.8p : 

• We introduce in (|1.2p an assumption of distributed fractional derivative (see |16| , 

a)- 

2 Fractional relaxation equation of order p 

A first probabilistic expression of the solution tp^{t),t > to equation p.2p can be 
found by considering that the latter coincides with the fractional equation satisfied 
by the survival probability (i.e. the probability of no events up to time t) of a frac- 
tional Poisson process of order i^ G (0, 1) . Let J\fu{t),t > 0, denote the process with 
probabilities p'j^ (t) solving the following recursive differential equation 



df 
with initial conditions 



'^'^^^^ -^Xipl^p'^,_,), fc>0,t>0, (2.1) 



^^(0) = ^ I l-\ (2.2) 



and p1]^(i) = 0. The process Mu has been studied in a series of papers (for example 
in [12], [m and [T], in the homogeneous case, and in [30], in the non- homogeneous 
case) and its distribution has been expressed in analytic forms in terms of derivatives 
of Mittag-Leffler function or as generalized Mittag-Leffler (GML) functions 

°° (7) 7J 
K^('^ = T. jir {aj+l3y «>/5,7eC, i?e(a),i?e(/3),i?e(7)>0, (2.3) 

where (7) ■ = 7(7-1- l)...{'j + j — 1) (for j ~ 1, 2, ..., and 7 7^ 0) and (7)0 = 1 (see [3]). 
Moreover in [T] a probabilistic expression of the process has been given, as composition 
of a standard Poisson process N with a random time argument 71- , independent of N. 
The following equality in distribution was proved to hold: 

Kit) ■=■ N{TAt)), (2.4) 

where Tu{t) possesses transition density qu{y, t) coinciding with the folded solution to 
the fractional diffusion equation 

" " i>0, yeR, i;(2/,0) = <5(y),«t(y,0) = (2.5) 



i.e. with 



-(»■') = < r^.'^o"-"- P-«) 



Alternatively, it has been also proved in [53| and in [TO] that q^{y,t) solves the fol- 
lowing equation 

|^ = -|> ^>0. 9(2/, 0)^%), (2.7) 

where, in this case, j/ > 0. In any case we can write 

Pl{t) = Pr{AC(t) = k}^-^J e~^yq,iy,t)dy, 



so that we immediately have, in view of p.ip for fc = 0, that 

Mt)^Poit) - Pr{AC(i) = 0} = / e-^yq,iy,t)dy. (2. 



Therefore, in view of (|2.4p . the fractional relaxation -01, can be expressed as composi- 
tion of the standard relaxation with the random time T^: 

2.1 Exponential boundary crossing probabilities of Brownian 
motion 



As a consequence of (j2.8p a second probabilistic form of the solution in terms of 
boundary crossing probabilities is obtained in the following result. 

Theorem 2.1 Let U be a random boundary exponentially distributed (with parameter 
A > 0), then the crossing probability of U by the independent random process T^{t) 
with transition density q^{y,t), i.e. 

Mt)^P^-{%{t)<U}, (2.9) 

satisfies the fractional relaxation equation il.2\) . with initial condition "0^(0) = 1. 
Proof We consider now the analytic expression of the folded solution qi,{y,t) to 
problem ()2.5p , in terms of the Wright function 

OO j 

W-/3(^) - E ,,r(a,+/3) ' " > -1, /? > 0, - e K, 

which reads 

q,{y, t) = 2v{y, t) = ^M/_,,i_, (-|;) , y, t > 

(see, for example, [Hj). Therefore we can rewrite (j2.9|) as follows 

?A,(i) = ¥T{Tu{t)<U) (2.10) 

= \ e-^y Pi- {%{t)<y}dy 



il?-'l>-A- 




hi:---—^^: 


)dz 


i?..i(-Ar), 





by the well-known formula of the Laplace transform of the Wright function (see |25| , 
formula (1.165), p. 39). The last expression in (|2.10p coincides with the solution to 
equation (|1.2p given in (|1.4p . ■ 



The previous results can be particularly relevant in the special case where v = 1/2, 
since the random process 71/ reduces to a reflecting Brownian motion: indeed in this 
case the equation ()2.5p governing the process coincides with the heat equation and 
51/2(2/, i) becomes the Gaussian with variance 2t, folded with respect to the origin. 
Therefore the fractional relaxation equation of order 1/2 is solved by 

1 r+°° 2 

^,/^{t) = ^= / e-''ye-'-dy^¥v{\B{t)\<U} (2.11) 



TTt Jo 

= Pr i sup B{s) < U 

[0<s<t 

The previous expression can be checked directly by applying ()1.4p : 

V'i/2(i) = Si/2,i(-AVt) (2.12) 

= [by the duplication property of the Gamma] 

- (-2AV^)^'r(i + i) 
^^ r(j + i)v^ 



1 r -. -1 ^ (-2xv^ty 



e ^z 5^. 



^-^0 - j! 



which gives (|2.1ip . after a change of variable. 

Also for u = 1/2", n > 1, the solution can be expressed in terms of boundary 
crossing probability of known processes. Indeed the random process 71/ coincides 
in this case with the (n — l)-times iterated reflecting Brownian motion defined as 
7n-i(t) — |i?i(|i32(---(|^n(t)|)---)l)li where Bj{t) are independent Brownian motions 
with variance 2f, for any j. The transition density qi/2" (y, t) of /„-i is given by 

2 ^2 ^2 

+00 r+°° e~^ e~^ ^ — ^ 



qi/2"iy,t) ^ ... ~l^--- — -^^dsi...dsn^i, 2/,t>0, 



which coincides with the folded solution to the following fractional diffusion equation 

rjl/2" o2 

■q^=q^^ yeR,t>0, g(y,0) = J(y), (2.13) 

(see [H], for n = 1 and [22], for n > 1). Therefore, in this case, the solution to 
the fractional relaxation equation can be expressed as the crossing probability of an 
exponential boundary by an iterated reflecting Brownian motion, i.e. 

/.+00 

^-1/2" it) = ^In-i (0) = / e-^^gi/2« (y, t)dy = Pr {/„_i (t) < U} . 
Jo 

For other rational values of the fractional order i/, such as. for example ly = 1/3, the 
solution can be still represented as boundary crossing probability, but of less known 
processes. 

For ly ^ 1/3 the random process % in (|2.9p reduces to the process A{t), introduced 
and studied in [22] , whose transition function is given by 

9i/3(y,0= Yy^*(^), y,t>o (2.14) 



where 



Ai{w) = - 



cos I aw + — — I da, w G 



(2.15) 



is the Airy function. By exploiting the relationship between (|2.15|) and the modified 
Bessel function 



/.H = E 



(f) 



2k+i' 



k=0 



k\T{k + l^+l)' 



w e 



I.e. 



Ai{w) 



1-1/3 



2Vw^ 



Ii/s. 



2Vw3 



w >Q 



we can rewrite the transition density (|2.14p of the process A{t),t > as 



91/3(2/, i) = \/^ 



^-i/^iVA)"^'/'rvi) 



y,t>Q. 



Therefore, in this case, the fractional relaxation can be written as 

/•+00 

V'i/3(t) = i>{A{t)) = / e-^^gi/3(y, t)dy = Pr {A{t) < U} . 
Jo 

It can be worth comparing the asymptotic behavior of the different crossing prob- 
abilities introduced so far. By using the well-known integral representation of the 
Mittag-LefHer function 



EuA-ctl 



t^-P 



^u-p -rt "t" sin(/?7r) -|- csin((/3 - v)tt) 
j,2v _j_ 2r'^ccos(;/7r) + c^ 



dr. 



we get the following asymptotic behavior of the solution ip^: 



Mt) 



1- 



\f 



At'-'r(i-i/) ' 



0<t << 1 

i — > 00 



(2.16) 



(2.17) 



Therefore the boundary crossing probability of Brownian motion exhibits a power 
decay, for t — > cxd, of exponent 1/2, instead of the usual exponential decay of the 
standard relaxation ip. For the n-tli times iterated Brownian motion the exponent 
1/2" of t is smaller than 1/2 and decreases as n becomes larger. This is intuitively 
explained by the fact that the number of subordinations increases in the definition 
of the process /„: this strays the fractional relaxation more and more away from the 
standard (exponential) behavior, as n increases, and makes the tail of the relaxation 
more and more heavy. 

For the process A{t) the crossing probability possesses a power decay, for t — > 00, 
with exponent 1/3 which is between the Brownian case and the iterated one (for any 
n> 1). 



2.2 



Exponential boundary crossing probabilities of more gen- 
eral processes 

We now present some extensions of the previous results, obtained by considering the 
crossing probabilities of different kinds of processes. This corresponds to substituting 
the random process T^{t) in p.9p with some other process, linked to the Brownian 
motion by various relationships, such as the elastic Brownian motion, the Bessel 
process (or its square), the first passage time through a level i by a standard Brownian 
motion or its sojourn time on the positive half line. 



We start from the latter, which, being a nondecreasing Levy process, can be con- 
sidered as a subordinator. Let F^ (t) = meas{s < t : B{t) > 0} be the sojourn time 
on the positive half line of a standard Brownian motion B, then its density g^(s, t) is 
given by 

ds 

g+(s,t) =Pr{r+ eds} = , 0<s<t. (2.18) 

ny/s{t- s) 

Theorem 2.2 Let U be a random boundary exponentially distributed, with parameter 
A > 0. Then the crossing probability of U by the random process T^ (t) with transition 
density g+(s,t), is given by 

^+{t) = ^(r+(t)) - Pr {r+(t) < [/} = e-^'/'-Io (^y) (2.19) 

and \2. 1 9]) solves the following second-order differential equation 

Proof We write the crossing probability as 

^+(i) = fe-^^ /; (2.21) 

Jo n^ys(t~s) 

= [formula 3.383.1, p.365 [5]] 
= iJ^iQ;l;-Ai 

where iFi (a,j;x) denotes the confluent hypergeometric function defined as 

^ a{a + l)...{a + j - 1) x^ 

iFi{a;r,x) = 1+ V— ^ — — ry-T — — 77-7, 

^ 7(7+l)-(7+J - 1) r- 

for cc, a € C and 7 G C\Z(^. 

By applying the relationship with the Besscl functions (see formula 9.215.2, p. 1086 
[9]) and, after some computations, we get the final form (|2.19p . As far as the equation 
satisfied by (|2.19p is concerned, we recall that /o(Aa;) coincides with the solution to 
the following equation 

^/o(Ax) + l^IoiXx) = X^IoiXx), (2.22) 

as can be easily checked. Therefore, by the transformation Iq (^) ~ e^*^^^p^{t), from 
equation (fZM|) wc get ^^ . since 

Alternatively wc can resort to to the form (|2.2ip and exploit the fact that the confiucnt 
hypergeometric function i.F'i (a;7;x) satisfies the following equation: 



d 
dx^ i i ' \ ' / ^2; 



^3Z2i^i + (7-^)-i^i=«i^i- (2-23) 



By taking into account that 

we get again (|2.20p . ■ 

The asymptotic behavior of "tp^it) can be deduced by considering that I^{x) ~ 
{x/2y/T{v + 1), as a; ^ 0, and that 

iFi (a; 7, x) ~ ^e-^^x"", Re{x) ^ -^ 
T{a) 

(see [H], p. 29), thus obtaining the following expressions 

L VAirt' 

The limiting behavior of ^^{t) is the same of a standard relaxation, for t — >■ 0, while 
coincides with that of '!/'i/2(^)i for i — ^ oo (up to multiplicative constants). 

Another process that can be considered instead of the random time Tu (t) in (|2.9p 
is the first passage time through a level i by a standard Brownian motion, denoted as 

T{t) =inf{.s >0 : B{s) = t} . 

Therefore, we are interested in the following crossing probability 



V'T(i) = Hnt)) 



/ e-^'qris, t)ds = Pr {T{t) < U} , (2.25) 

-'0 



where the density of T{t), i > is the well-known stable law of index 1/2, i.e. 

,.(M) = ^^, M>0. 

Therefore (|2.25p can be easily evaluated, since the Laplace transform of the first 
passage time is well-known: 

^^{t) = e-*^^. (2.26) 

Clearly ipxit) satisfies the standard relaxation equation, even if with a different con- 
stant: 

^=-V2A^T, ^^(0) = !. 
at 

We remark that time-changing the relaxation V' by the 1/2-stable subordinator T{t) 

produces again a standard relaxation, while performing the same operation by the 

inverse stable subordinator E{t) wc get the fractional relaxation -01 (as mentioned in 

the introduction). 

If we now consider n independent Brownian motions Bj, j = 1, ...,n and construct 

by them the n-times subordinated process Ti(T2{...Tn{t)...)), t > 0, where Tj = 



inf {s > : Bj{s) = t} , j = 1, ...,n, then its crossing probability can be evaluated as 
follows: 



Pr{Ti(r2(...r„(t)...)) 


<c/} 


(2.27) 


/ e-^^ / dz,. 
Jo \Jo 


r-^., te-V... ,„_ 




-io "^""-^^2^- 


V2^z3 V27rs3 y 


r+oo /• + 00 

/ dzi... / dZn- 
Jo Jo 


ie-*'/2-i z„_ie---"-i/2^' 




' V2t^? ^27^4 


Jo V27rs3 


.+00 te-tV23i 


/•+°° z„_ie-^"-/2^" ^_,„ 


^^dz 1 


Jo y/2Trzf 


io v/27rz3 




.+00 ie-*'/2-i 


/■+°°z„_2e-^"-/2.„_,^_ 




3„_iV2V2A^_, 


Jo VStt^i 


^" V^"""-! 




g-A2Tr2l-2Trt 







Again the probability i/ii^ satisfies (for any n) the standard relaxation equation with 
the constant A 2^ 2^^ 2^ and displays an asymptotic behavior similar to the standard 
relaxation, despite the complicated construction via the n-times subordination. 

We analyze now the crossing probability of an exponential boundary ?7 by a squared 
Bessel process. Let us denote by i?^(i) = (i?-y(i)) , < > the square of a 7-Bessel 
process, starting at zero. It is well known that, for 7 = n, this process can be expressed 
as 

n 

Rl{t)=Y,B^{t), i>0, 

where Bj(t),j = 1, ...n, are independent Brownian motion in M". Moreover the density 
of i?2 can be written as 



^^^^'^^ = (2t)ir(i) ' ^'*>" 

(see, for example, [7]), which is a more tractable form (for our aims) than that of R-y. 
Thus the crossing probability of this process can be easily evaluated as follows: 

i>^{t) = Pr{i?2(t)<t/} (2.28) 

_;^5 8 2 ^e 2t 

(2t)*r(i) 
1 



(2At + l)^ 
which satisfies the following first-order differential equation 

1^- = 2A^'^-' ^-^°) = ^- 

In this case, the behavior of ip^yit), for increasing (but still finite) values of i, 
can be represented as ij-yit) ~ (fc/i)^/^ (for some constant k and for < 7 < 2) 
and thus it coincides with the one described as "algebraic decay" and displayed by 
relaxation processes in complex material (see, for example, |27]). On the contrary, for 
the other fractional relaxations, this is true only in the limit, for t — >■ 00. Indeed the 



function (|2.28p coincides with the so-cahed Nutting law, which is commonly used to 
fit the experimental data for the materials featuring non-standard (i.e. non-Debye) 
relaxation (see [IH] and the references therein). 

As we have seen, the generalizations analyzed so far in this section are not linked 
to fractional equations; on the other hand, in the following case, we consider crossing 
probabilities governed again by fractional equations. Let B"{t), i > be the so called 
elastic Brownian motion with absorbing rate a > (see ^U\ and [2]), defined as 

Bf(i) = { '''^'^'' ^^JJ" , (2.29) 

where T^ is a random time with distribution 

Py{T^>t\Bt} = e-°'^^°'*\ a > 0, (2.30) 

Bt = (t{B{s),s < t} is the natural filtration and L(0,i) = lime4,o^"T-eas{s < t : \B{t)\ < e} 
is the local time in the origin of i?. It is well known that its distribution can be ex- 
pressed as 

q^J{s,t) = 2e"' we-""" -^=dw + qc,{t)5{s), s, t > (2.31) 

J s v27ri3 

where 5{s) is the Dirac's Delta function with pole in the origin and 

1 ,..2 

q^{t) = 1 - Vv{B''^{t) > 0} = 1 - 26^^ / '-j=dw 

JaJt v27r 



'a\/t v27r 

is the probability that the process is absorbed by the barrier in zero up to time t. 
Thus we define the crossing probability of an exponential boundary U by the process 
Bf! as 



C'(i)=Pr{i?S'(i)<t/}= / 



e" C(s,Orfs- (2.32) 



Theorem 2.3 Let U he a random boundary exponentially distributed, with parameter 
A > 0. Then the crossing probability of U by the random process B^{t) with transition 
density q^{s,t), is given, for any X ^ a, by 



rj{t)^PT{B^j{t)<u}^i ^ 



A — a 
while, for a = A, it coincides with 



--l-^)--(-f) 



<(i) = PriBfit) < [/} = 1 - XV2iE.,. (^-^) 



, (2.33) 



(2.34) 



The crossing probability ipai^) satisfies, for any a, A > 0, the following fractional 
differential equation 



dt^''^ V2 dii/2^" 2 ^ ^"^ V2^' 



>"+^:^^"=^(i-^^')-:7fe' ^S'(o) = i. (2.35) 



10 



Proof We take the Laplace transform of (|2.32p . which reads, for any a, A > 0: 



Jq Jo 



(2.36) 



2 / e-''*dt I e 
10 Jo 

/•OO 



-Xs+as 



ds 



we 



e 2t 

-.dw 



^/2^ 



1-2 / e-^'+^dt '-^dw 



V Jo 



2 / e 
/o 

2 






V 



2a 



1 /■+°° 1 

e 






2T] + a Jo 
2 



ri(2ri — a'^) 



a{y/27i~ a) 

(^ + a){y/2^ + X) ^ v{2tj - a^) 

_ aXij^'^ + V2ar]^i +2 

Wc can check that p.36p coincides with the Laplace transform of p.33[) . for a ^ X, 
as follows: 

/•OO 

Jo 

1 r / \ i / \ \ il foo 

e-'^Hidt 



1 


V 


1 


( "V (^ ^ V 


77 


A-a^ 

3= 


1 

7? 


A 1^ 

A - an ^'^ 


Y a y / A y" 




1 


A 1 


■ \/2^ \/2^ ■ 
V2^ + a ^/2^+A 




T] 


A — a ry 


' 



which easily gives (|2.36p . As a further check of (|2.33p . it is easy to see that, for a = 
(in the case of no absorption) it reduces to ipi{t) = i?i/2,i(— A\/f), since in this case 

B"'(i) = |B(i)|,i>0. 

For a = X the Laplace transform p.36p becomes 



°°e-Vi'(0^t^-^'^"+^^^"^+2 







(V2^ + A)2 



(2.37) 



By comparing (|2.37p with the formula holding for the Laplace transform of the GML 
function defined in (|2.3p (see [II], p. 47), i.e. 



C{t^-'El^{u:t0);r^] 



0/35-7 



1 

(where i?e(/5) > 0, -Re(7) > 0, Re{S) > and 7] > \uj\ ^-f^) ), we easily obtain 



rx\t) = i 



V2 2-2 
11 



V2)' 



(2.38) 



(2.39) 



which can be also rewritten as (|2.34p . 

By taking the Laplace transform of equation (|2.33p and considering the well-known 
expression for the Laplace transform of the Caputo derivative, i.e. 



dt" 



-r,t 






i{t)dt 



(2.40) 



m-l ,r 

= ^^C{u-,i^}-Y,if-<'-'—u[t) 



t=o 



we get 






A 



aX A 
~2\ 



^/2 

Ar(i) 

y/2Trq 



ri"^C{rI\n] 



-^^-^rJ^iQ) (2.41) 



By taking account the initial condition "0^(0) = 1, the solution of (|2.4ip coincides 
with (P^Sj) . ■ 

In order to study the asymptotics of the solution ^p^{t), for a 7^ A, we use the 
integral expansion for the Mittag-LefHer function (|2.16p . so that we get 



rJit) 



1 



A 1 



A-: 



+ 00 



z-i/2e-- 



V2 



A 

V2 



'■-Vt 



f. + 'tVi, 



L Vt ^ 2 
Therefore the limiting behavior of the crossing probability reads 



dz. 



rJit) 



XV2t 

V2 
Q\/7ri ' 



< t « 1 

t — ^ 00 



(2.42) 



(2.43) 



where the first line is obtained from (|2.42p by the following calculations: 

xVt '■+°° 



cw 



1 



V2t:Jo 
V27r 



z~^/^e~'dz 



[by the reflection formula of Gamma function] 
XV2i 



Thus, in this case, the crossing probability maintains a limiting behavior similar 
to the previous ones for i — > 0, but drastically different for i — > 00 (see ()2.17p ). 
In the last case instead of tending to zero, it tends to one: this can be intuitively 
explained by considering that the absorbing effect is stronger as t increases and, in 
the limit, the process B"^' will be absorbed with probability one. This effect is directly 
correlated with the absorbing rate a. Thus it is evident from (|2.43p that ip^ looses 
the usual property of complete monotonicity that characterizes the standard and also 
the fractional relaxations (see, for example, |16j). 

In the case a = X we must apply the integral expansion of GML functions (see 

a) 



E'sh^tn 



,l^P fOO 



2711 



g -rt^uk-P 



^inp 



12 



dr, 



(2.44) 



(for k = 2, ly = 1/2, j3 = 3/2 and c = A/a/2) so that formula (|2.39p can be developed 
as 



<(i) 



= 1 



1 



A 1 



V2 2TrJo { 



Af.^2 
2 / 



V2 



A 1 



"Z 2 



ii + ^r 



dz. 



V2 



dr 



/2tTT JO 

Therefore, also for a = A, the asymptotic behavior is given exactly by ()2.43p . 



Remark 2.1 An interesting link can be found between the crossing probabilities 
ip^{t) and ipi/2{t)'- for A 7^ a, the first one can be rewritten, in view of (|2.33p and 
(ETH), as 

A — a I ' ' 



rjit) 



(2.45) 



where ip'^i^it) and '0i/2(^) indicate the crossing probability Pr{|i?(t)| < [/} of an 
exponential boundary U of parameter a and A, respectively, by a Brownian motion. 
Thus the following identity is also verified for the corresponding differential equations: 



d}^ 

dty^ 



r: = 



A 



A — a 
A 



dl/2 



dt 



dl/2 



TTI^V^W- 17171^1/2(0 



A — a 
by applying Theorem 2.1, for j/ = 1/2 



^/2 






V'?/2(i) - 7^^1/2(0 



2.3 Crossing probabilities of a Gamma distributed boundary 

Wc extend the previous results by considering the crossing probabilities of a random 
boundary, distributed with different laws, instead of the exponential one. In particular 
we choose its natural generalization, i.e. the Gamma distribution. Thus we are 
considering the following probability, which extends formula (|1.8p 



/•OO 

V'i(i)=Pr{|B(i)|<G}= / [l-Faiy)] 
^ Jo 



o-vVit 



-dy, 



(2.46) 



where G is a Gamma r.v. with parameters A, A; > and Fq denotes its cumulative 
distribution function. For our convenience, we write the latter as follows: 



Fciy) 



A* 



nk) Jo 



'-^^z^-idz=(^^) 



r(fc) fr'.jKj + ky 



■E 



(2.47) 



Theorem 2.4 Let G be a random boundary distributed as a Gamma with parameters 
A, fc > 0. Then the crossing probability of G by a standard Brownian motion is given 
by 

^'Ht) - Pr{|B(t)| < G} = 1 - {XVifE'! . , l(-A^/t), (2.48) 



which satisfies the following fractional relaxation equation 



(2.49) 
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with initial condition i/'llO) = 1, for k>l, and the additional conditions 






0, r = l,..., 



t=o 



t=o 



, /or any orfrf fc > 1 



0, r = 1, ..., — — 1, /o?' any even k > 2. 



Proof Wc can rewrite ()2.46p as 



V'Kt) = 



k °o 



Et 



i-xyy 






= -tfV4t 



-dy 



(-iyA-'+^ p .+, 



r(fc)V^^ j!a + ^-) Jo 



E^-t^HV/ y^-^'^e-^/^'rfy 



^ 00 



{-i)H2xViy+^ ^ f 3 fc 1 



r(^)V^,t^ j!(j + fc) 



222 

2 ^ (-ip(Ax/tp'^ r(j + fc-) 



rW^t^ jiU + k) r(| + |) 
(Av^)'= ^ r (j + fc) i-xVty 



(2.50) 



If wc now assume that fc is an integer, we can recognize in ()2.50p the GML function 
p.3p . so that we get p.48p . As a further check, it is easy to ascertain that, in the 
special case fc = 1 (where the r.v. G reduces to the exponential r.v. [/), the crossing 
probability ip't given in (|2.48p coincides with the fractional relaxation ipi in (|2.12p : 



ipi{t) = l-XVt.Eisi-XVt) 

= 1+E ^'n^K ^E.,{-XVt)=^PL{t). 



(2.51) 



In order to derive equation (|2.49p we resort to the Laplace transform of (j2.48p which 
reads: 

i^+xr-x'^ 



£ 



{^|;,} 



^(V^+A)* 



(2.52) 



by applying again formula (|2.38p . for 7= -1 + 1, /3= ^ and S = k. Wc now rewrite 
(|2.52p as follows 



E 

3=0 



Xk~3 



i]iclxjA;7]\ -ri'^ ^ 



X^ 



(2.53) 



By simplifying this expression, we can recognize the Laplace transform of equation 
(|2.49p . We can check that the initial conditions are satisfied, by using the series 
expression of i?^_^(— ci''), and considering that for t = 0, E^j^{—X^/i) = l/r(/3): thus 
we get 

(Ax/t)^ 



mt) 



= 1- 



t=0 



r(| + i) 



= 1. 



t=0 
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For the other conditions, we can apply the following formula of the r-th order deriva- 
tives of a GML function (see formula (1.9.6), p. 46 of [TT]): 






z''-'Elp{\z") 



^P-r-l^P 



a,P — r 



(Az"), 



A e C, r e N, 



so that we get 



cf 






(2.54) 



(2.55) 



By recalling formula (|2.40p , we notice that the Laplace form (|2.53p holds if the deriva- 
tives of order r of f/'i vanishes for r ~ 1, ... [IJ if fc > 1 is odd and for r = 1, 



■•2 



.1 
if fc > 2 is even; this is verified by (|2.55p . 

Finally we check that equation (|2.49p becomes, for fc = 1, the fractional relaxation 

equation -^^iH) — ~\ibi(t). ■ 

dt2 2^ ' 2^ ' 

Remark 2.2 By comparing (|2.48l) with the results in [3], we can deduce that the 
crossing probability il}\ (t) can be written in terms of the fractional Poisson process 

2 

of order i^ = -i , as 



xp\{t) = Vv{Tk >t} = Pr {aAi (t) < fc} 



(2.56) 



where T^ = inf -^ t > : A/i (i) = fc > is the waiting probability of the fc-th event. On 
the other hand we can prove that the following relationship holds between the crossing 
probabilities given in (|2.46l) for a Gamma boundary of parameters (A, fc) and (A, fc — 1) 
(respectively denoted as ip\{t) and -0j^^^(t)): 



di/2 



dt 



-L/'^ 2 L 2 2 



(2.57) 



Indeed we can evaluate the fractional derivative of order 1/2 of ip\, by considering 



di/2 



dt 



\^ 



^^l^ = -7^(fc^§ 



(j + fc - l)!(-A)^" 



,=0 J!r (i + f ) Jo 



{t-sY 



-^ds 



(2.58) 

Z 'Z ■ z 

By applying to p.58p the following recursive formula for GML function proved in ([3]) 
^"^™„.+.(--^) + 2:"+'K'(„+i).+,(-x) = x"i?™„-V.(-2;), n, m > 0, z > 0, X > 0,^ 
for TO = n = fc, a: = — A\/f, i/ ~ 1/2, z = 1/2, we can rewrite 



dti 



= -t-^ IxHiEt-'.i^xVt) - x^+h'^E^ ..,hxV~t) 

U'T"'" 2 2 ' 2 "'"-^ 

= -A'=i*-^i?^l , (-AVf) + A(l - i'lit)), 

2'2 + 2 2 

which gives p.57p . The latter could be alternatively obtained by considering that 
pl^'it) = PrJAA. (f) = fc} = V4(t) - i^l^'it) 



(2.59) 
(2.60) 
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satisfies (|2.ip with v ~ 1/2 and taking into account (|2.56p . 

The asymptotic behavior of the crossing probability -0^ for small t can be deduced 

2 

by the series expression of the GML function 

^tA-^n^^^^^y 0<t«l, (2.61) 

so that we get 

The same result can be obtained by resorting to the Laplace transform and to the 
Tauberian theory, which permits to infer (formally) the asymptotic behavior of a 
function /(t), for i — ^ oo and t — > 0^, from the limiting behavior of its Laplace 
transform £ {/; r]} for r/ -^ 0+ and t] ^ oo, respectively (see also [H], for details). To 
this aim, we rewrite (|2.52p as 

which, for ry — > cx3, can be approximated as follows 

cUhv} = --^+oiv-^-') (2.64) 

SO that we get again f|2.62p . For t — > oo, it is worth writing f|2.52p as 

so that we get -0^ it) ~ ^, — . Thus the limiting behavior of -0^ can be summed up as 
follows: 

V'i(t)^<! ^W^' , (2.65) 

which, of course, coincides with ()2.17p for fc = 1 and i' = 1/2. We can deduce that, 
while for small t passing from an exponential boundary to a Gamma-distributed one 
makes a relevant difference, for large t this effect fades away. Indeed the rate of the 
decreasing to zero for t — ^ oo of the crossing probability is exactly the same for any 
fc > 1. 

Analogously, we can generalize the results of Theorem 2.3, by considering the 
crossing probability of a Gamma distributed boundary by the elastic Brownian motion 
defined in (jT^ . 

Theorem 2.5 Let G be a random boundary distributed as a Gamma with parame- 
ters A, fc > 0, then the crossing probability of G by the random process B'^{t) with 
transition density q'^^{s,t) (given in V2.31\ ). for any A, a > 0, is equal to 

(2.66) 
16 



which, in the particular case a = A, reduces to 



rUt) = ^^{Bf(t) <G]=l- (^) E'^l^^ (-^) . (2.67) 

Proof By following some steps similar to those of Theorem. 2.3, we can write the 
Laplace transform of i/jf^it) as follows 



-CK^;^}-/ e-^'dt [l-FG{s)]qt!{s,t)ds (2.68) 

Jo Jo 

pCC I' + OO 1 

2/ [l~FGis)]e°"ds e-("+^^)"'du> + - - 

Jo Js V 

2 2a 



2ry-a2 ^{2ri ~ a"^) 

2 2A2 ^/fc + j-1 



(V2^ + a)V2^ ^/2^''+'(V2^ + a)^t^ V J 
rj{2ri — oi^) 



+ 



V2^(V2^+a)(V2^+A)'= ?7(V2^ + a) 
1 V2A^ 

We can invert (|2.68p by applying again (|2.38p : 

- 1 ,/2A^£f ^ ^''^' .f 

^'^''\(V2^ + a)(y2^+A)'='' 



^ -(;!)/(- ^)-^^^-^-^(-^)--H... (-^)- 



V2; ^r(i + i)^(fc-i)!j!r(i + Mi) 



X 



/ / s 1 — i fc-1 I J , 

X / (i- s)2 ^s 2 +2c?5^ 

which, after some simplifications, coincides with (|2.66p . For a = A, we can rewrite 
the latter as follows: 

-^) 

^^^^ i^,iz(fc-i)K"i-0!r(f + 1 + 1) 

^ (\^t\^^ [-^T ^(k + m-l-l\ 

VV2J ^^^r(f + 1 + 1)^1 m-l ) 
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= [by the identity proved in 0, p. 10] 
= 1- 






m=0 V 2 ' 2 

As a final check, we can ascertain that, for fc = 1, formulae p.66|) and p.67p reduce 
to the corresponding expressions given for the exponential case in (|2.33p and (|2.34p . 
respectively: indeed (|2.66p can be rewritten, for fc = 1, as 



<«w = 1+E -^ E 



oo ( aVt\ oo / /--. I oc ( ^Vi\ 



^/2; t;^v V2; fr^ r(| + i) 

oo f_Ay5^ fc 

^. r f ^ + 1) f- 



a\^\ ,v^l^ V^'J V^/^o^A' 



which coincides with (|2.33l) . Formula (|2.67l) immediately reduces to the expression 
tirm . for fc = 1. 

Finally, putting a ~ and substituting A/\/2 with A, formula (|2.66p coincides 
with the corresponding crossing probability ()2.48p . which has been obtained in the 
case of a free Brownian motion (with no absorption). ■ 

The asymptotic behavior of V'Jc'a' f°^' small t, can be derived from (|2.67p . by ap- 
plying again formula ()2.6ip . Alternatively we can use the Laplace transform (|2.68p . 
which can be approximated as follows, for rj ^ oo 

V 22 7/2+-^ 

In both ways, we get the first line of the following formula: 



r,Ut)^{ Mv^J r(i+i)' o<t«i ^^^^^ 

The second line of the previous expression has been obtained from (|2.68p . which can 
be rewritten as 



'' VV [E -=0 ( ■)(2'7)^+^A-^ +aE -.0 ( ■)(2r;)*A-^ 
1 V2 






For fc = 1, formula (|2.69p coincides with (|2.43p . as was expected. We finally note 
that, also in this case, as for the Brownian motion, the leading term in the expression 
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obtained for i — >■ oo does not depend on k and thus, for large values of t, considering 
an exponential or a Gamma distributed boundary does not entail any consequence. 

The fractional equations satisfied by the crossing probabilities obtained above can 
be derived by properly rewriting the Laplace transform in p.68p . as the following 
theorem shows. 



Theorem 2.6 The crossing probability V'^'q given in 



satisfies, for any A, a > 0, 



the following fractional equation 



" 'k\ fV2y di 



Ck 



sc) (fit^j't..^m m s*'-^(-<)-.. 



(2.70) 
where Ck = 1 for k odd and Cfe = for k even. The initial conditions are fp'^'ai^) ~ ^' 
for any k > 1 and 






t=0 



t=0 



= 0, r = 1,. ..,——, for odd k>l (2.71) 



0, r = l, ..., — — 1, for even k > 1. 



Proof We rewrite (|2.68p as follows: 



so that we 


get 










(V2\ 


<.^^+^ 


a 

V2 


'1 

.V 


'V~^:„ 


- 


1 



T^Efnf^Vfe.*-.*-' 



a 

72 



j=i 



(2.72) 



where we have denoted ipk a " '^ i ^fc'aj "H \ ^^r brevity. From the Laplace transform 
(|2.72p . by taking into account (|2.40p and the initial conditions (|2.7ip . we can obtain 
equation ()2.70p with c^ = 1. For the initial conditions ()2.7ip we use an argument 
similar to that of Theorem 2.4, with the only additional care that, in the case of even 

k 

k, the highest order derivative, i.e. -^V'fe'ai does not vanish in t = 0, as can be 
ascertained by applying (|2.54p to ()2.66p : indeed we get 



at 2 



1 \ fe oo 

^)5 



-^\ t'.FM 



V2j 



t^E'l 



xVt\ 

V2j 



V2 



Therefore formula (|2.72p . for even k, must be modified as follows 

fe 



m(4)' 



3=0 



el ^i + 



^kcV 



2-1-1 2_1 

2 "'" 2 ^ 712 2 



V2' 



V'fea'7' 



1 2 2 



7f 5*-<" 



t=0-i 



a sr-^ fk\ j v2 



A 



7el i i-1 



a 

71 



1 



^r. 
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so that we get (|2.70p . with Cfc = 0. 

As a further check, it is easy to see that, for fc = 1, the latter reduces to equation 

3 Fractional relaxation equation of distributed or- 
der 

Wc consider now an extension of the fractional relaxation equation ()1.2p obtained 
by adding the hypothesis that the fractional order v is not a constant but a random 
variable with distribution n{v). Thus we will study the distributed order fractional 
relaxation equation defined as 

J -^n{iy)diy ^ ~X^, i > 0, (3.1) 

where, by assumption, 

n(z/)>0, / n{v)diy = l, i^ € (0, 1] , (3.2) 



subject to the initial condition ip{0) — 1. As a special case, for n{i') ~ Sliy — v) and a 
particular value of F S (0, 1) , equation p.ip reduces to (|1.2p . 

We adopt here the following particular form for the density of the fractional order 
v. 

n{v) = ni5{v — vi) + n25{v — V2), Q <vi <V2 <l, (3.3) 

for 71,1,712 > and such that ?t,i + 77,2 = 1 (conditions p.2p are trivially fulfilled). 
The density p.3p has been already used by [T7] and [5], in the analysis of the so- 
called double-order time-fractional diffusion equation, and corresponds to the case 
of a subdiffusion with retardation. Moreover, it was applied in j3] in the context 
of recursive equations of fractional order, where the equation governing the Poisson 
process has been extended by introducing two fractional time derivatives. 
Under assumption p.3|) . equation (|3.ip becomes 

d"^ d"^ 

and the corresponding solution iIjvi,v2 coincides with the so-called double- order frac- 
tional relaxation studied by [11]. They provide for ipvi,w2 a-n integral expression and 
some asymptotic representations. We present here an analytic form of the fundamen- 
tal solution to p.4p in terms of GML functions as well as a probabilistic representation 
in terms of crossing probabilities, in line with the results of the previous sections. 



Theorem 3.1 The solution to equation \3Mj with the initial condition ■(/;(0) = 1 can 
he written as follows: 

Proof By taking the Laplace transform of p.4p we get 

ni77''i£{V'j.i,i.2;'7} -7?"' +"277'''-C{V'i^i,i.2;^}-'7''' = -A£ {^/'i.i,,.^;^} : (3-6) 
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whose solution can be written as 



riiX + niri"^ + n2if^) 
1 A 1 1 

1 A 1 



?; A + n2r]"2 1 + "i''"^,, 

' ' X-\-7i2r]'^ 



■E 



?/ ?/ A + 112^]"^ ^ \ A + n2rf'^ 



n\r\ 



■E 



'? "2 — 



"2 



^v\r — \ 



' n2 



r+1 



By applying formula (j2.38p . wc easily get (J33|). As a check we can see that p.Sp 
reduces to (|2.10p . for ni = 0, n2 = 1, j^2 = '^i since equation p.4p becomes, in this 
case, the fractional relaxation equation (|1.2p . ■ 



Despite the apparent similarity of p.Sp with (|2.66p . they are deeply different: while 
for "01^1, 1/2 the sum is extended to the third (upper) parameter of the GML function, 
this is not the case for f/'I-'a- This is also reflected in the asymptotic behavior of 
the fractional relaxation of distributed order, which does not deviate from the usual 
relaxation behavior (unlike i^k^)- We can study the limit directly from (j3.5p . by 
applying formula (|2.44p . as follows 



A ■^-^ fTli\ 1 






"2 ^ V "-2/ 27ri Jo 



g-^t^i.ir-1 



(3.7) 



^— z7ri/2^i'?T{i/2^t'l)^ 



giixVi+iT^iv-i—vxjr 



"2 



r+1 



r+1 



Thus, for t — > 0, we get 

n2f^\ "2 / 27ri Jo 



e-^y^-^ir-i 



(3.8) 



- -i TT 1/2 — i TT { ^'2 — 1^ 1 ) ^ 



p'i7ri/2+'i7T(i/2 —J^i)f 



\ "2 



\J-1'2 °° 

ri2 ^ 
r— 



r+1 



^^'^■^ ^i-KUo 



r+1 



711 i'^ 



sin(— 7r(i/ir — 1^2^ — ^2)) 



712 / TT 

= [by the reflection property of the Gamma function] 



r {vir - V2r - V2) 



r=0 



n2 



At" 



1 



r(l +7/27' + 1^2 - I'l?') 



772 r (1 + V2) 



+ o(r=). 
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while, for t — >■ cx), we analogously have that 



A 



ni 



r—0 ^ 



^ — z7ri/2 — Z7r(;>2~^i )^ 



1 

27ri 



^-w^.,r-l 



(3.9) 



^i7n^2-+-'i'^{i^2—i^i)r 



r+1 



Y^ / ni ^^ sm(7rz^ir) 

r=0 

oo 



r (z^ir) 



A pi-nU2 



r+1 



r=0 



At'^i 



1 



ni 



1 



r(i-j^ir) Ai^^i r(i-i^i; 



+ o(t-'^0- 



The previous expressions coincides with formula (4.16) of J16j . which has been ob- 
tained in a different way, directly from the Laplace transform of 'ipui,u2- 

We present now a probabilistic form of the solution ip^-^^^^, which is in line with 
the analysis carried out so far, in terms of crossing probability of a random boundary 
by a stochastic process, that will be denoted, in this case, by 7I'i.i/2(t),i > 0. To this 
aim we will compare equation p.4p with the equation governing the probabilities pk 
of the distributed order fractional Poisson process Mui.u2{'t):'t- > studied in [1], i.e. 



^d'^Pk 



dt" 



n{v)dv = -X{pk-Pk-i), fc>0, p_i(t)=0 



(3.10) 



Indeed p.ip can be considered a special case of p.lOp for fc = and, if we add the 
assumption p.3p . we get p.4p . Thus we can use the results proved in [3] and write 
that 



V'.L.aW =Po(i) = Pr{AA,,,,,(t) = 0} = ¥r{N{%,^,,{t)) = 0} 



(3.11) 



where N is the standard Poisson process (with intensity A) and Tui,u2 is a random 
process (independent from N) with density 

qut,u2{yit)^ni I p^^{t- s;y)q^^{y,s)ds + n2 p^^{t - s;y)q^^{y, s)ds. (3.12) 
Jo Jo 

In p.l2p Pj^.(-;z) denotes the density of a stable random variable X^. of index i^j e 

(0, 1] , for j ~ 1, 2, with parameters equal /3 = 1, /i = and a = (?^j|y| cos ^^) ' 
and q^ , for j = 1, 2, was defined in ()2.6p . Another form of the density qv^^u2 is given 



by the following series expression 










<l:^l,^2iy,t) 








= 


oo _, , 
r=0 ^ 


"2|y|\ 
XV"^ ) 






we get 


oo 

XtU2 L^ ^\ 
r=0 


f nilyl 
\ At'^i 


From 


(|3.11| 








"01^1, i 


,M= r '~ 


^Vn ( 

'i'^l,'^2\ 



(3.13) 



W- 



1/1,1 — 1-2^ I/l 



y^-V2,l-vir-U2 



ni\y\ 
At'^i 

n2\y 



Xt''^ 



{y,t)dy^Vi{%,,u2{t) <U}. 



(3.14) 
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It is also proved in [1] that the transition density qui,v2 coincides with the folded 
solution 

of the following fractional diffusion equation 



"i^^+n2^^j =^, yeM,t>0, ni,n2>0, (3.16) 

for < 1^1 < i^2 < 1, with initial conditions 

v{y, 0) = 6{y), for < z^i < i^2 < 1 .„ -, -j-. 

wMy^t)\t=o = for i < i/i < ^2 < 1 ■ ^ ' 

In alternative to p. 16^ - ^.171) it can be proved (as we will see below in a special case) 
that qui,u2 solves also the other equation 

Ql^ly Ql'2y Qy 

m-^^+n2-^^ = -^, y,t>0, ni,n2>0, v{y,0)=S{y), (3.18) 

which is the distributed order analogue of (|2.7p . In order to get a more explicit 
expression of the density qui,u2^ we consider the special, but relevant, case where 
1^1 = ^ and 1^2 = 1- 

Theorem 3.2 The solution to the fractional relaxation equation 

"i^^+ri2^ = -AV^, t>0, (3.19) 

with the initial condition "0(0) = 1, can be expressed as follows: 

^i,i(t)=Pr{ri,i(t)<[/}, (3.20) 

where U is an exponential r.v. with parameter A and the transition density of Ti iit),t > 
0, is given by 

2 2 

gi,i(2/,0- ^^ r- ^, i>0, 0<y<-, (3.21) 

and satisfies the fractional equation 

^^Sft+"^l = -S' '^(^.0)=%). (3.22) 

Proof It has been proved in [4] that for 1^2 = 1 and vi = v ^ (0, 1) the density 
(|3.12p . can be expressed as 

quAy,t) = nir{f^{-;y)){t) +n2p^{t;y), (3.23) 

where /'' is the Riemann-Liouville fractional integral of order v and p^^ denotes a stable 

law of index 1/ and parameters equal to /3 = 1, /i = 'T-2|y|, f = (^T-ilyl cos 2^) . If 
we put moreover i> = 1/2, we can recognize in pi the Levy distribution, so that the 
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density p.23p becomes 

= — L / {t- sy'^fi{s]y)ds + n2piit]y) 
V^r Jo ' ^ 



27r 



(i-,s)- 



g i(.B-n2V) 



zds 



nin2y e ■!i(*-"2a) 



2tt 







2^ v/(i-n22;)3^{°<^<4} 



[t — n2y — z) 2 — ___q(2 + 



[by the identity (3.8) of [2T] 



nie 4(-'-'»2!/) 7117122/ e ''('-"ay) 



^7r(t - nay) 27r ^{t-n2yf 



\^<y<^y 



which coincides with p.2ip . In order to show that the latter satisfies the fractional 
relaxation equation p.22p . we evaluate its Laplace transform, which reads: 






(3.24) 




\J{t-n2yY 



dt 



-r,t 



^ g 4(t-,i2H) 

y ^/{t-n2yY 



rdt 



nin2y 
2V^ JO 

-Tin-2y-^niy 



-r)Z-7)7l2l/ ^ ' ^„ 



= (772+77l?7-l/2)e-("^''+""'''')^ 

In p.24p we have applied the well-known formula of the Laplace transform of the 
first-passage time of a Brownian motion. It is easy to check that 



-^^£{gi^i(y,.);?/}rfy 



?i2 + niT] ^''^ 

7727? -I- 7^1771/2 -f A' 



which is equal to the Laplace transform of i}}iJi,v2i for ^1 = 1/2 and !^2 = 1 (given in 
Theorem 2.6 of jl]), thus proving result p.20p . If we now take the Fourier transform 
of ([3:241) we get 



J". 



{£{9i,i;77};/3}- H e^^y c[q^_^,{y,■);,^] dy (3.25) 



772 + "1?? 



-1/2 



^iPy ^-{n2V+niri^''^)y ^ 



772 + nit] ^ " 
7727/ + 7717/1/2 + if3^ 

which coincides with the solution to equation p.22p converted, via Laplace-Fourier 
transform, into 

(77177I/2 + n2T])C [qi.iiy, •); '7} - {ni'n^^''^ + n2)5{y) = ~-Q-^ {9ii(y' ■)'''} 
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and 

From p.25p it is evident that p.2ip is well-defined and integrates to one, since for 
/3 = we get 1/ry. ■ 

Remark 3.1 If we consider the two opposite special cases 7^2 = and rii — 0, the 
trajectories of the process 71 1 can be considered as "interpolation" between those of 
a free reflecting Brownian motion and the straight line y = t/n2. Indeed in the first 
case the density p.2ip becomes 

2 2 
nie it 



while in the second we can write p.23p as qi i{y,t) = n2pi(t]y) = n2(5(i — n2y), 
since in this case cr = 0. It is evident from (|3.2ip that the trajectories of 71 1, for any 
7X1,712 > are forced under the line y — t/n2 and this is reflected in the asymptotic 
behavior of the crossing probability i/ii j, which can be deduced from ()3.8p and p.9p 
and summed up as follows: 

(^ X\/izt 

By comparing p.26p with ()2.17p we can conclude that i/ji i displays the same limiting 
behavior of 7/'i(t) ~ Pr{|_B(t)| < [/}, for t — > 00. On the contrary, for t ^ 0, it 
behaves as the standard relaxation (up to a constant) and thus tends to one much 
faster than ipi i. We recall that similar limiting features were exhibited by the crossing 
probability ■0+ of the Brownian sojourn time process (see (|2.24p ). 

For the reader's convenience we summarize the limiting behavior of the crossing 
probabilities analyzed in the previous sections in the following tables: 



Table 1: Limiting behavior for f — > 



iPit) ~ 1 - Ai 
Mt) ^ 1 ^*' 



r(i+z.) 
0+(i) ^ 1 - ^ 
V'^(f)~l- V2Ai 
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Table 2: Limiting behavior for i — >■ oo 

yt^W — AtT(l-iy) 

"1 



V'i.iW 



\sfni 
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